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Abstract 

We present a measurement of the Cabibbo-Kobayashi-Maskawa CP-violating phase 7 with a 
Dalitz analysis of neutral D-meson decays to the KgTr~7r~^ final state from — > D^*^K^ de- 
cays, using a sample of 347 million BB events collected by the BABAR detector. We measure 

7 = (92 lb 41 lb 11 ± 12)°, where the first error is statistical, the second is the experimental system- 

(*) 

atic uncertainty and the third reflects the Dalitz model uncertainty. For the ratios between the 
magnitudes of amplitudes A{B^ D^*^^K~) and A{B^ D^*^^K^) we obtain the one-standard 
deviation intervals [0,0.14] and [0.02,0.20], respectively. All results presented here are preliminary. 
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1 INTRODUCTION 



The angle 7 of the unitarity triangle is the phase of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix [T| defined as 7 = 8i,ig[—V^jV*fj/V^^V*i^], which corresponds to the phase of the element 
V*f^, i.e. Vub = \Vub\e~^'^ , in the Wolfenstein parameterization |2j. Various methods have been 
proposed to extract 7 using D^K^ decays, all exploiting the interference between the 

color allowed B~ — > D^K~ {b ^ us (x Vcb) and the color suppressed B~ — > D^K^ (6 ucs oc 
Vub) transitions [Sj, when the and are reconstructed in a common final state |^ El E]. 
The symbol indicates either a or a meson. The extraction of 7 with these decays is 
theoretically clean because the main contributions to the amplitudes come from tree-level diagrams 
(see Fig. lU . 



U 




Figure 1: Diagrams contributing to i? — > D K decay. The left diagram proceeds via b ^ y,s 
transition, while the right diagram proceeds via b — > ucs transition and is color suppressed. 



Both BABAR 'W,^^ and Belle jTUj have reported on a measurement of 7 based on i? D^*^^K 
and B- D^K*' decays with a Dalitz analysis of D° K'^Tr-TT+, with L>°7r° and 

D*o ^ Z)°7 (BABAR only), and K*' K'^n' . In this paper we report on an update with 
B- bi*)^K- decays. 

Assuming no CP asymmetry in D° decays, the B^ i)^*^^K^, Z)°7r°,L)°7, 
Kg7r~TT~^ decay chain rate T^\m?i,m'^) can be written as [B] 

r^;\ml,ml) oc \ADT\^ + r^B^'\AD±\^ + 

2kr^*^ {cos(4*^ T 7) MADTAh±] + sin(4*^ T 7) MADTAh±]} , (1) 

where and are the squared invariant masses of the Kg'ir~ and Kgir^ combinations, re- 
spectively, and An^f = ADifn^jiri^), with Ad- {Ad+) the amplitude of the —>■ Kg7:~7r~^ 
(D^ Kg'K~^Tr~) decay. The value of the CP-odd phase 7 changes sign for B'^ and B~ in Eq. (pQ), 

leading to different rates in corresponding regions of the Dalitz plane, for B^ and B^ decays. 

(*) (*) ■ - ' 

We introduce here the CP (cartesian) parameters and defined respectively as the real 

and imaginary part of r^^e*'-'^s'^'''\ for which the constraint = xfj*^ holds. Here, is 

the magnitude of the ratio of the amplitudes A{B- D'^*^^K-) and A{B- L>(*)°i^-) and is 
their relative strong phase. As a consequence of parity and angular momentum conservation in the 
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Z)(*)° decay, the factor k in Eq. O takes the value +1 for D^K^ and D*°{D°tt^)K^ , 

and -1 for B^ D*^{D^j)K^ HI]. 

Once the decay amphtude Ad known, the Dahtz plot distributions for from B~ 
D'-*^°K- and B+ D^*)^K+ decays can be simultaneously fitted to T^*\m'^, m^) and T^^\m'^ , m"^) 
as given by Eq. (pQ), respectively. A maximum likelihood technique is used to measure the CP- 
violating parameters ■ From them, confidence regions for 7, and 6^^^ are obtained with 

a frequentist method. We extract y^'^ instead of 'J^S^b^ ,r^B^ because the distributions of the 

(*) (*) 

cartesian parameters are unbiased and Gaussian, while the distributions of 7,5^ don't have 

(*) 

these properties for small values of and low-statistics samples. 

2 THE BABAR DETECTOR AND DATASET 

The analysis is based on a sample of approximately 347 million BB pairs collected by the BABAR de- 
tector at the SLAG PEP-II e~^e~ asymmetric-energy storage ring. The BABAR detector is optimized 
for the asymmetric-energy beams at PEP-II and is described in 12.. We summarize briefly the 
components that are crucial to this analysis. Gharged-particle tracking is provided by a five-layer 
silicon vertex tracker (SVT) and a 40-layer drift chamber (DGH). In addition to providing precise 
space coordinates for tracking, the SVT and DGH also measure the specific ionization (dE/dx), 
which is used for particle identification of low-momentum charged particles. At higher momenta 
{p > 0.7 GeV/c) pions and kaons are identified by Gherenkov radiation detected in a ring-imaging 
device (DIRG). The typical separation between pions and kaons varies from 8<t at 2 GeV/c to 2.5a 
at 4 GeV/c. The position and energy of photons are measured with an electromagnetic calorimeter 
(EMG) consisting of 6580 thallium-doped Gsl crystals. These systems are mounted inside a 1.5 T 
solenoidal super-conducting magnet. 

3 EVENT SELECTION 

We reconstruct the B" b^*)^K- decays with ^ D°TT^,D°-f and D° K'^tt-tt+ 0. 

The Kg candidates are formed from oppositely charged pions with a reconstructed invariant mass 
within 9 MeV/c^ of the nominal Kg mass 113] • The two pions are constrained to originate from 
the same point. The — > i^^7r~7r''' candidates are selected by combining mass constrained Kg 
candidates with two oppositely charged pions having an invariant mass within 12 MeV/c^ of the 
nominal mass Jl3|. The vr'^ candidates from D*^ — > D^n^ are formed from pairs of photons 
with invariant mass in the range [115, 150] MeV/c^, and with photon energy greater than 30 MeV. 
Photon candidates from D*^ D^j are selected if their energy is greater than 100 MeV. D*^ 
D^TT^{D^^) candidates are required to have a D*^-D^ mass difference within 2.5 (10) MeV/c^ of its 
nominal value ^^1) corresponding to about two standard deviations. B~ — > D^*^^K~ candidates 
are formed by combining a D^*^^ candidate with a track identified as a kaon. 

We select B mesons by using the energy difference AE = E'^ — E* /2, and the beam-energy 

substituted mass, m-Es = {E*"^ /2 + pi • p-^Y /Ef — p\, where the subscripts i and B refer to 
the initial e^e~ system and the B candidate, respectively, and the asterisk denotes the center- 
of-mass (GM) frame. The resolution of AE ranges between 15 MeV and 18 MeV depending on 
the decay mode. The resolution of rn-Es is about 2.6 MeV/c^ for all the B decay modes. We 
define a selection region through the requirement —80 < AE < 120 MeV and rriES > 5.2 GeV/cP. 
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To suppress e+e" qq, q = u,d,s,c (continuum) events, we require |cos07^| < 0.8 where 6t 
is defined as the angle between the thrust axis of the B candidate and that of the rest of the 
event. Furthermore we define a Fisher discriminant J- that we use in a likehhood fit to separate 
continuum and BB events. It is defined as a Unear combination of four topological variables: 
Lq = J2iPiJ -^2 = J2iPi \ cos0*p, the absolute value of the cosine of the CM polar angle of the B 
candidate momentum, and | cos^tI- Here, p* and 9* are the CM momentum and the angle of the 
remaining tracks and clusters in the event, with respect to the B candidate thrust axis. If both 
B^ D*^{D^'K^)K^ and B~ — > D*^{D^'y)K~ candidates are selected in the same event, only the 
B^ D*^{D^iT^)K^ is kept. The cross-feed among the different samples is negligible except for 
B^ D*^{D^^)K~ , where the background from B~ D*^{D^7r^)K~ is about 5% of the signal 
yield. This contamination has a negligible effect on the measurement of the CP parameters. 

The reconstruction efficiencies are 15%, 7%, 9%, for the D^R-, B~ L>*°(£»°7r°)K- 

and B~ — > D*^{D^^)K~ decay modes, respectively. Fig. |21 shows the tjies distributions after all 
selection criteria plus a tighter requirement on AE, \AE\ < 30 MeV, are applied. The largest 
background contribution is from continuum events or BB decays where a fake or true is com- 
bined with a random track. Another source of background is given by those B~ D^*'^'^tt~ decays 
where the prompt pion is misidentified as kaon. These decays are separated from the signal using 
their different AE distribution. 

4 The L»o ^ A^^Tr 7r+ DECAY MODEL 

The KgTT~'TT~^ decay amplitude ^z)(m^,m^) is determined from an unbinned maximum- 

likelihood fit to the Dalitz plot distribution of a high-purity (97.7%) sample from 390328 D*^ 
D^ir^ decays reconstructed in 270 fb~^ of data, shown in Fig. El Our reference model to describe 
Aoirn'^, 771^) is based on Breit-Wigner (BW) parameterizations of a set of resonances, and is the 
same as used for our previously reported measurement of 7 on B- D^*)^K-, B- D^K*-, 
Z)° ^ K'^TT-Tr+ decays j3El. 

The decay amplitude in the reference model is expressed as a sum of two-body decay-matrix 
elements (subscript r) and a non-resonant (subscript NR) contribution, 

AD{ml,ml) = Erttre''^- Ar{rr?__,m\) + aNRe^'^^R ^ (2) 

where each term is parameterized with an amplitude a-r (onr) and a phase ('/'nr)- The function 
.4.r.(m^,m^) is the Lorentz-invariant expression for the matrix element of a D'^ meson decaying 
into i^g7r~7r~'' through an intermediate resonance r, parameterized as a function of position in the 
Dalitz plane. For r = p(770) and /3(1450) we use the functional form suggested in Ref. while 
the remaining resonances are parameterized by a spin-dependent relativistic BW distribution. The 
angular dependence of the BW terms is described with the helicity formalism as shown in |15j^. 
Mass and width values are taken from ^nj; with the exception of i^o(1430)"'" taken from ^Hj. The 
model consists of 13 resonances leading to 16 two-body decay amplitudes and phases (see Table 
plus the non-resonant contribution, and accounts for efficiency variations across the Dalitz plane and 
the small background contribution. All the resonances considered in this model are well established 
except for the two scalar tttt resonances, a and a', whose masses and widths are obtained from our 
sample |17j . Their addition to the model is motivated by an improvement in the description of the 
data. 

^The label A and B should be swapped in Eq. (6) of |15| . 
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Figure 2: Distributions of mss for (a) D^K', (b) B- D*°{D°Tr^)K- , and (c) ^ 

D*^{D^-))K~ . The curves superimposed represent the overall fit projections (solid black lines), 
the continuum contribution (dotted red lines), and the sum of all background components (dashed 
blue lines). 



The possible absence of the a and a' resonances is considered in the evaluation of the systematic 
errors. In this respect, the K-matrix formalism ^Sj provides a direct way of imposing the unitarity 
constraint that is not guaranteed in the case of the BW model and is suited to the study of broad 
and overlapping resonances in multi-channel decays. We use the K-matrix method to parameterize 
the TTTT S-wave states, avoiding the need to introduce the two a scalars. A description of this 
alternative model can be found in [S]. 
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Component 






r It 11 action {/o) 


K*{892)~ 


-1.223 ± 0.011 


1.3461 ± 0.0096 


58.1 


i^^(1430)" 


-1.698 ± 0.022 


-0.576 ±0.024 


6.7 


i^2*(1430)" 


-0.834 ±0.021 


0.931 ±0.022 


3.6 


K*{14:W) 


-0.248 ± 0.038 


-0.108 ± 0.031 


0.1 


ii^ (1680) 


—1.285 ± 0.014 


0.205 ± 0.013 


0.6 


K*{892)+ 


0.0997 ± 0.0036 


-0.1271 ± 0.0034 


0.5 


i^^(1430)+ 


-0.027 ± 0.016 


-0.076 ± 0.017 


0.0 


i^|(1430)+ 


0.019 ± 0.017 


0.177 ± 0.018 


0.1 


p(770) 


1 





21.6 


u;(782) 


-0.02194 ± 0.00099 


0.03942 ± 0.00066 


0.7 


/2(1270) 


-0.699 ± 0.018 


0.387 ± 0.018 


2.1 


p(1450) 


0.253 ± 0.038 


0.036 ± 0.055 


0.1 


Non-resonant 


—0.99 ± 0.19 


3.82 ± 0.13 


8.5 


/o(980) 


0.4465 ± 0.0057 


0.2572 ± 0.0081 


6.4 


/o(1370) 


0.95 ±0.11 


-1.619 ±0.011 


2.0 


a 


1.28 ±0.02 


0.273 ± 0.024 


7.6 


a' 


0.290 ±0.010 


-0.0655 ± 0.0098 


0.9 



Table 1: Complex amplitudes a^e*'^'' and fit fractions of the different components {Kstt~, Kstt'^, 
and 7r"'"7r~ resonances) obtained from the fit of the Kstt~'7t'^ Dalitz distribution from — > 

D^n~^ events. Errors are statistical only. Masses and widths of all resonances are taken from jl3) 
with the exception of -R'q (1430)"*" taken from ^B]. The fit fraction is defined for the resonance terms 
as the integral of af.\Ar{m'^,m'^)\'^ over the Dalitz plane divided by the integral of |w4i:)(m?_, m^)p. 
The sum of fit fractions is 119.5%. A value different from 100% is a consequence of the interference 
among the amplitudes. 



5 CP ANALYSIS 

We simultaneously fit the D^*^^K^ samples using an unbinned extended maximum-likelihood 

fit to extract the CP-violating parameters x)^ and along with the signal and background yields. 
The fit uses tti-es, AE, T ^ and . The likelihood for candidate j is obtained by summing the 
product of the event yield the probability density functions (PDF's) for the kinematic and 
event shape variables Vc^ and the Dalitz distributions "pDaiitz^ ^^^^ signal and background 
components c. The likelihood function is 

/: = exp[-^iV,WE^^^^(C)^?^''^(^.) ' (3) 

\ c / j c 

where = {mes, AE, f]j = {ml,ml)j, and 3(1) = 7^i,c(mEs)7'2,c(AS)7'3,c(^). The 
background components in the fit are continuum, BB and B^ D^tt~ (for B~ — > D^K~) or 
B~ D*^TT~ (for B~ D*^K~). For signal events, V^''^^^^^ (ff) is given by r^\f]) multiphed by 

the efficiency variations estimated using simulated signal events, where r^^(?7) is given by Eq. 

The rriEs and AE distributions for signal events are described by Gaussian functions; the Fisher 
distribution is parameterized with two Gaussian functions with different widths for the left and 
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m? {GeV^/c") ml,. (GeV^/c") 

Figure 3: (a) The Kg-K^n^ Dalitz distribution from D*^ D^tt^ events, and projections 

on (b) = m^o^+, (c) = mj^o^-, and (d) m^+^_. i^^vr+vr^ from D*^ D^tt^ 

events are also included. The curves are the reference model fit projections. 



right parts of the curve (bifurcated Gaussian). Their parameters, along with most of the parameters 
describing the background distributions, are determined from a combined fit to the B~ D^*^^7t~ 
high-statistics control samples. 

5.1 Description of the background probability density functions 

The continuum background in the mss distribution is described by a threshold function jl9| whose 
free parameter is determined from the B~ D^*^'^Tr~ control samples. The continuum AE 
distribution is described by a first order polynomial whose slope is extracted from the control 
samples. The shape of the background m^s distribution in generic BB decays is taken from 
simulated events and uses a threshold function to describe the combinatorial component plus a 
bifurcated Gaussian shape to parameterize the peaking contribution. The fraction of the peaking 
contribution is extracted directly from the fit to the data. The AE distribution for BB background 
is taken from simulation and parameterized with the sum of a second order polynomial and a 
Gaussian function that takes into account the increase of combinatorial feed-down background at 
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negative AE' values. The m-Es distribution of — > D^*^'^Tr~ is the same as the signal, while 
the AE shape is parameterized with the same Gaussian function as the signal with an additional 
shift arising from the wrong mass assignment to the prompt track, computed event by event as 
a function of the prompt track momentum in laboratory frame and the CM boost. The Fisher 
PDF for continuum background is determined from the m-Es sideband region of the control sample 
events and is parameterized with the sum of two Gaussian functions. The Fisher PDF for BB 
events and B^ D^*^^tt^ background is taken to be the same as that for the signal, consistent 
with the simulation. 

Background events arising from continuum and BB where the candidate is real can mimic 
either the 6 — > c or the b ^ u signal component, depending on whether the candidate is 
combined with a negatively or positively-charged kaon. We take this effect into account in the 
likelihood function with two parameters, the fraction fj^o of background events with a real 
and the fraction R of background events with a real associated with a negatively-charged kaon 
(same charge correlation as the b ^ c signal component). These fractions have been estimated 
separately for continuum and BB backgrounds from simulated events. As a check of the reliability 
of these estimates, the fraction f^o for all background events (mixture of continuum and BB) has 
been measured on data from the invariant mass distribution of after removing the requirement 
on the mass and using events satisfying rriEs < 5.272 GeV/c^. The measured value is consistent 
with the fraction found on simulated events. The fractions fjjo and R for continuum and BB 
background are reported in Table [21 

The shape of the Dalitz plot distribution of the continuum and BB background is parameterized 
by a third-order polynomial function in (m?_,m^) for the combinatorial component (fake neutral D 
mesons), and as signal or shapes for real neutral D mesons. The combinatorial distributions 
are taken from simulated events. The shapes for events in the invariant mass and rriEs sidebands 
on data and simulated events are found to be consistent. The fraction of background originating 
from signal B~ — > D^*^'^K~ where the D^*^^ meson is combined with a combinatorial (either 
opposite- or same-charged) kaon from the other B meson is found to be negligible. 



fraction 


B- D^R'- 




B- D*^{D°j)K- 


fj^o (continuum) 
R (continuum) 
/do (BB) 
R (BB) 


0.022 ±0.010 
0.164 ±0.018 
0.026 ± 0.008 
0.64 ±0.15 


0.336 ± 0.038 
0.170 ±0.052 
0.130 ±0.041 
0.5 ±0.5 


0.163 ±0.016 
0.099 ± 0.031 
0.152 ±0.024 
0.943 ± 0.039 



Table 2: fractions f^o and R, as described in the text, from simulated continuum and BB 
background events. 

5.2 CP parameters 

The signal yields measured with the CP fit on the sample of 347 million BB events are N{B^ — > 
Z>0/s:=F) = 398±23, iV(5=F ^ Zj*°(l)°7r°)KT) = 97±13, NiB"^ D*°{D°-y)K'f) = 93±12, and are 

consistent with expectations based on measured branching fractions and efficiencies estimated from 

(*) (*) 

Monte Carlo simulation. The results for the CP-violating parameters , are summarized 
in Table El The only non-zero statistical correlations involving the CP parameters are for the 
pairs (x_,y_), (x+,y+), (a;!.,y*), and (x!j_,y!j_), which amount to -1%, 1%, -17%, and -14%, 
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respectively. The Dalitz plot distributions for the events selected with rngs > 5.272 GeV/c^ are 
shown in Fig. |1] separately for and candidates. Fig. |H1 shows the one- and two-standard 
deviation confidence-level contours (including statistical and systematic uncertainties) in the x^*^ — 
y^*^ planes for all the reconstructed modes, and separately for B^ and B^ . The separation of the 
{x^-\y- '^) and {x^_^\y^^) confidence contours in these planes is an indication of direct CP violation. 



CP parameter 


5=F ^ £)(*)0X=F 


X- 


0.041 ± 0.059 ± 0.018 ± 0.011 


y- 


0.056 ± 0.071 ± 0.007 ± 0.023 


x+ 


-0.072 ± 0.056 ± 0.014 ± 0.029 


y+ 


-0.033 ± 0.066 ± 0.007 ± 0.018 


x*_ 


-0.106 ± 0.091 ± 0.020 ± 0.009 


y*- 


-0.019 ± 0.096 ± 0.022 ± 0.016 


xX 


0.084 ± 0.088 ± 0.015 ± 0.018 


y+ 


0.096 ± 0.111 ± 0.032 ± 0.017 



Table 3: CP-violating parameters , obtained from the CP fit to the B^ — > samples. 
The first error is statistical, the second is experimental systematic uncertainty and the third is the 
systematic uncertainty associated with the Dalitz model. 

5.3 Systematic error associated with the Dahtz model 

The largest single contribution to the systematic uncertainties in the CP parameters comes from 
the choice of the Dalitz model used to describe the — > Kg7r~TT~^ decay amplitude. The D'^ 
sample used to determine the reference model introduced in Secl^is fitted with a set of alternative 
models where the resonances are described with different parameterizations or removed: 

1) vrvr S-wave: the reference model uses two wide BW scalar amplitudes {a and a'). Alternatively 
we use a K-matrix model Pl with pole masses and coupling constants fixed by fits to scattering 
data [201 • See also Sec. 01 

2) TTTT P-wave: the mass and the width of the Gounaris-Sakurai BW describing the p{770) are 
changed within their quoted uncertainty ,13 . 

3) TTTT and Ktt D-waves: alternative to the helicity formalism used in the reference model, for 
/2(1270) and -fC|(1430) we use the formalism derived from Zemach tensors [21j. The difference 
is very small for P-waves but is larger for D-waves. 

4) Kit S-wave: the mass and width of the BW describing i('*(1430) are taken from E791 |16j . 
Alternatively, we have floated them in our flavor tagged sample obtaining consistent 
values. As an additional model we use an adaptation of the LASS parameterization |22l with 
parameters taken from the fit to our D*~^ D^tt^ data sample. 

5) Ktt P-wave: it is dominated by the K*{892) in both Cabibbo allowed and doubly Cabibbo 
suppressed amplitude. The mass and the width of this resonance, taken from PDG jl3| in 
the reference model, are changed to the values found by keeping them floating in the fit to 
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Figure 4: The L>° ^ i^"°7r-7r+ Dalitz distributions for D^K^ (a,b), 5=F ^ ^*0(Z)07rO)is:=F 

(c,d), and — D*^{D^^)K^ (e,f), separately for (a,c,e) and i?"*" (b,d,f). The requirements 
mES > 5.272 GeV/c^ and lAE'l < 30 MeV have been apphed to reduce the background contami- 
nation. 
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Figure 5: Contours at 39.3% (dark) and 86.5% (light) confidence level (corresponding to two- 
dimensional one- and two-standard deviaton regions), including statistical and systematic uncer- 
tainties, for the {x^\ y^^) parameters for (thick and solid lines) and (thin and dotted lines) 
decays. 



the flavour-tagged D sample. The resulting values are consistent with what is found in 
B J/'ifK-K decays selected in BABAR data. 

6) Blatt-Weisskopf penetration factors: the effect from the Blatt-Weisskopf penetration factors 
has been evaluated using an alternative model that doesn't include them [2'6\ . 

7) Running width of BW: a model with BW's of fixed width is used. 

8) K2{14:30), K*{1680), i^*(1410) and p{1450): these resonances are removed from the reference 
model. 

We have generated a sample of B^ D^K^ and B^ D*^K^ signal events that is one hundred 
times larger than the measured signal yields in data. The Dalitz plot distribution of is generated 
according to the reference model and to CP parameters consistent with the values found in data. 
The CP parameters are extracted by fitting the generated Dalitz plot distributions using a PDF 
equal to the reference model (model 0) or to one of the eight alternative models {model 1, 2,...,8). 
We take as the systematic uncertainty of {xzf,yzf) — similarly for (xj^,y^) — associated with the 
i''^ alternative model the difference between the CP parameters fitted using the alternative model 
(x^,yy and the reference model {x^,y9^): Ax*p = — x!^, Ay?j- = y*p — y^- As total systematic 
uncertainty associated with the Dalitz model we consider the sum square of contributions from 

the alternative models: Axzp = J J^^^i Ax'^J^ , Ay^ = J'}2^=i Ay'^J^ . The dominant contributions 
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to the overall Dalitz model uncertainty arise from the models 1), 4), and 7). The systematic 
uncertainties associated with the Dalitz model are summarized in Table 0J 



Soiivrp 


x_ 


y— 






X* 


U— 


X*, 


y+ 




0.002 


0.004 


0.003 


0.004 


0.011 


0.012 


0.008 


0.008 




002 


000 


n 000 


n 000 


n 002 


n 003 


002 


n 01 6 


Fraction of right sign D^^s 


0.008 


0.002 


0.002 


0.002 


0.005 


0.005 


0.001 


0.022 


Efficiencv in the Dalitz plot 


0.014 


0.000 


0.013 


0.001 


0.001 


0.002 


0.000 


0.001 


Background Dalitz shape 


0.006 


0.003 


0.001 


0.004 


0.012 


0.015 


0.009 


0.009 


Dalitz amplitudes and phases 


0.004 


0.004 


0.004 


0.004 


0.008 


0.008 


0.008 


0.008 


B D*^K- cross-feed 


0.000 


0.000 


0.000 


0.000 


0.004 


0.001 


0.004 


0.004 


CP violation in Dtt and BB bkg 


0.000 


0.000 


0.000 


0.000 


0.005 


0.002 


0.002 


0.005 


Total experimental 


0.018 


0.007 


0.014 


0.007 


0.020 


0.022 


0.015 


0.032 


Dalitz model 


0.011 


0.023 


0.029 


0.018 


0.009 


0.016 


0.018 


0.017 


Total 


0.021 


0.024 


0.032 


0.019 


0.021 


0.027 


0.023 


0.036 



Table 4: Summary of the main contributions to the systematic error on the CP parameters x=p, y^, 
x^, and y^. 

5.4 Experimental systematic errors 

The main experimental systematic errors are listed in Table 0] Uncertainties due to the m-ES) 
AE, and PDF parameters for signal and background extracted from the combined fit to the 
B~ — > D^*^^'K~ control samples (fixed in the reference CP fit) are estimated from the statistical 
differences on and when the former set of parameters is also floated in the CP fit. Other 
rriES) AE, and J- parameters fixed in the CP fit are changed by one standard deviation. The 
uncertainties associated to the knowledge of the fraction of background events with a real 
and the Dalitz distribution of background events are evaluated from the differences on the CP 
parameters when the estimates obtained from simulated events are replaced by the estimates using 
sideband data. The systematic uncertainty on the fraction of events where a true is associated 
with a negatively-charged kaon is obtained from the variation of the CP parameters when the 
is randomly associated either to a negatively- or positively-charged kaon (absence of charge 
correlation). The effect due to reconstruction efficiency variations of the signal across the Dalitz 
plane has been estimated assuming a perfectly uniform efficiency. The statistical errors in the 
Dalitz amplitudes and phases from the fit to the tagged sample have been propagated to the 
xfj*^ and parameters performing a simultaneous CP and Dalitz fit to the B^ D^*^^K~ and 
D*^ —>■ D^TT^ data. The effect of the remaining cross-feed of B^ D*^{D^'k^)K~ events into the 
B~ D*^{D^'y)K~ sample (5% of the signal yield) has been evaluated by including an additional 
background component with 'P^^''*^(r/) identical to that of B~ — > D*^{D^n^)K^ signal events. 
Finally, possible CP-violating effects in the background have been evaluated by setting the CP 
parameters of the B~ — > D^*^^'k~ background component to the values obtained from a CP fit to 
the B~ D^*^^'ir~ control samples, and by fioating an independent set of CP parameters for the 
other BB background. 

The following sources of uncertainty are found to be negligible: the assumption of perfect mass 
resolution for the Dalitz plot variables (m^,m^), the presence of combinatorial background from 
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signal events where the prompt kaon is replaced by a combinatorial track, and the assumption that 
the shape of the continuum or BB background does not change when the is fake or real. 

6 INTERPRETATION 

A frequentist (Neyman) procedure ^1 IS| identical to that used in our previous measurements [HI 
IH] has been adopted to interpret the measurement of the CP parameters {x^\y^^) reported in 
table 01 in terms of confidence regions on p = {'^ b , ^ b if*^-, ^*b) ■ Using a large number of pseudo- 
experiments with probability density functions and parameters as obtained from the fit to the 
data but with many different values of the CP parameters, we construct a multivariate Gaussian 
parameterization of the PDF of ( function of p which takes into account the statistical 

and systematic correlations. For a given p, the five-dimensional confidence level C = 1 — a is 
calculated by integrating over all points in the fit parameter space closer (larger PDF) to p than 
the fitted data values. The one- (two-) standard deviation region of the CP parameters is defined 
as the set of p values for which a is smaller than 3.7% (45.1%). Figure IHl shows the two-dimensional 
projections onto the rB—j and r^—j planes, including statistical and systematic uncertainties. The 
figure shows that this Dalitz analysis has a two- fold ambiguity, (7,5^^) ^ (7 + 180°, + 180°), 
as expected from Eq. (^. From the one-dimensional projections we obtain for the weak phase 

7 = (92 ± 41 ± 11 ± 12)°, and for the strong phase differences 6b = (118 ± 63 ± 19 ± 36)° and 
6^ = (— 62ib59±18±10)°. No constraints on the phases are achieved at two standard deviation level 
and beyond. Similarly, for the magnitude of the ratio of decay amplitudes vb and we obtain the 
one (two) standard deviations constraints vb < 0.140 (r^ < 0.195) and 0.017 < r|j < 0.203 (r|j < 
0.279). All these results are obtained considering the statistical correlations mentioned in Sec. 15.21 
while the experimental and Dalitz model systematic uncertainties are taken uncorrelated. We 
have verified that accounting for experimental systematic correlations within a given measurement 
(xzp,yzp) or (x?j-, yj!j-), or assuming the experimental and Dalitz model systematic uncertainties 
between (x^,y^) and {x^,y^) fully correlated, has a negligible effect on the results. 

7 CONCLUSIONS 

We have presented a preliminary updated measurement of the CP parameters (xzp, r/zp) and {x^, y^) 
with S=F ^ D°7r^,D°j, D° K'^tt-7t+ decays based on a data sample of 347 

million BB pairs, that supersedes the previous one based on about 227 million BB pairs i^j. The 
current analysis reduces the experimental systematic uncertainty and improves the procedure to 
estimate the error associated with the Dalitz model of the decay. 

Despite the improved measurement of (x^ , ) , the uncertainty on 7 has increased with 
respect to our previous measurement [H], moving from 7 = (70 ± 3llj^oiii)° to 7 = (92 it 41 it 
11 lb 12)°. Since the uncertainty on 7 scales roughly as l/r^^\ this change is explained by noticing 

that the new {x^\y^^) measurement favors values of smaller than our previous analysis and 
significantly smaller than the latest Belle results [TUj . 
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Figure 6: Projections in the (a) — 7 and (b) — 7 planes of the five-dimensional one- (dark) 
and two- (light) standard deviation regions. 
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